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SECTION I 
INTRODUCTION 


The original version of the SINGER program was developed under 
contract F29601-73-C-0022. Since the completion of that phase of the 
development, significant experience has been gained working with the 
program. Avariety of problems has been studied, including those with 
steel elements and those with reinforced concrete elements. Much of 
the detailed understanding acquired has been documented for future 
reference, and code modifications have been introduced to correct 
errors or to improve computational accuracy. The objective of this 
volume is to document these areas of study for future reference. 

The organization of the documentation is in the following form: 

1. Additional user assistance for the preparation of input 

data and for the interpretation of results; 

2. Additional program details not available in previous 

documents; 


3. New subroutine documentation; 


4. Listing of the program code in its current form. 


The objective of this section is to describe alterations and 
clarifications in the preparation of problem data for input and in the 
explanation of some output details consistent with the current status 
of the program. The input and output features of the original SINGER 
code were both documented in the User's Guide [1]. A broader 
understanding of the program has provided the basis for expanding this 
documentation for the benefit of the user. Input data preparation is 
discussed in Addenda to the User's Guide and Error Controls. Some 
output characteristics are discussed in Minimization and Convergence 


Problems and Strain Discontinuity. 


ADDENDA TO THE USER'S GUIDE 

Although the basic structure of the original User's Guide is still 
valid, some points require correction or clarification. In addition, 
a few modifications need to be noted. These items are listed below; 
the page numbers shown refer to the original User's Guide [1]. 

1. Corrections: 


a. 


defines the characteristic of the drop-elastic 


SECTION II 
ADDITIONAL USER ASSISTANCE 


(p. 23), footnote no. 3: The parameter k, which 


unloading curve for concrete, must be greater than 
zero or less than or equal to 1 to be an acceptable 
value. 

tp. 34), footnote no. 4: If the user inputs the 


number of bars and the bar size number, subroutine 
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BARS will calculate the area, bar size diameter, 


and perimeter. 


c. (p. 41), footnote no. 1: If no initial conditions 
are specified and IANAL = 0, a static analysis is 
performed according to specified load increments. If 


no initial conditions are specified and IANAL = 1, a 


static analysis is skipped and a dynamic analysis is 
performed. 
d. (p. 42), card 3/R, columns 1 - 5: Format I5 should 
be format 14, 1X. 
e. (p. 43), footnote no. 3; paragraph 3: For vertical 
F | Slopes, i.e., t;,, - t; = 0 (where t represents time), 
the forcing function value is set to F; in subroutine 
TABL. 
2. Clarifications: 
a. (p. 18), footnote no. 6: The maximum relative error 
for determining the converged state of minimization 
t | (SERR) cannot be less than vEPS, where currently EPS = 
1.E - 14; otherwise an error message is printed. If 
this value is input as zero, the program sets its 
value to VEPS. 
b. (p. 22-23), Material Data Block: User input material 
functions must be defined by seven points, each with a 
corresponding stress and strain value. None of these 


strain values can be in descending order. The concrete 


material must have negative stress and strain values, 
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but the crushing strength should be a positive number. 

The stress value at point 5 is used to indicate 
a default function for confined concrete for stress and 
strain modification due to confinement effects. If 
this stress value is zero, the stress and strain values 
at point 5 and the strain at point 6 are computed; if 
this stress is nonzero, a user input function is 
indicated and no additional computations are needed. 
(Although this computation affects the material 
function, it is not done in the input sequence until 
after the stirrup data are known.) 

Both unconfined and confined concrete material 
functions should be input for a problem using reinforced 
concrete elements. If a single concrete material is 
desired for the element, both functions can be input 
with identical values. 

The steel stress-strain values input by the 
user should be positive. 

The strain corresponding to the peak stress of 
the user input function for concrete is the reference 
point for determining the unloading characteristics; 
unloading prior to reaching this reference strain value 
is a straight line parallel to the initial segment 
slope (path a b in figure 1); unloading after passing 


this strain value follows the drop-elastic path 


determined by the parameter k, where 0 < k 2 1.0 


I eg 


FIG. 


PLASTIC OFFSET STRAIN 


1 


CONCRETE MATERIAL PARAMETER DEFINITIONS 
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(path c d in figure 1). The drop-elastic parameter 
k is used to determine the length of the vertical stress 
segment as shown in figure 1. 

The array S containing nine elements is used to 
keep track of stress history data in energy density 
computations for concrete. These are defined in 
subroutine CRET. The first eight elements, S(1) thru S(8) 
are called the stress history data, while the ninth 
element, S(9), is referred to as the plastic strain 
data. (All values are stored in the DATA array for 
each element of status 3.) The eighth element, S(8), 
is called the effective plastic offset strain value, 
and this quantity is defined specifically in figure 1. 
(A similar set of nine values is defined for steel 
material in REIN.) 

Several input values of material constants 


are checked to see if they fall within a specified 


range. 
0.= [Poisson's Ratio] <0.50 
0.8 [Shear modulus] <15.x10° psi 
0.£ [Material density] <0.50 1b/in° 
0.< [Unloading constant] £1.0 
2.5x10°psi < Unconfined concrete strength 5 Sg. x10°psi 


“Confined concrete strength 


33.x10°psi £ [steel yield stress] £75.x10°psi 


area Re a 


Cc. 


(p. 28-39), Element Data Block: Each reinforced concrete 
element requires the element parameter card and the concrete 
data card. It ts advisable to have at least one 
longitudinal reinforcement group since no check is made 

to see if there is none, and some DO loops are executed 
over the number of groups. The lateral reinforcement 

card is optional. (Stirrups do not enter directly into 

the element energy computations.) An unreinforced section 
can be simulated by introducing a very small nonzero 

steel reinforcement group area. 

Each steel wide flange element requires only the 
wide flange reinforcement card; and a leaf spring element 
(if used) requires only the leaf spring flexibility 
card. 

The lateral reinforcement parameter AY ip. 35) 
refers to the cross sectional area of the stirrup at 
a single location, i.e. the cross sectional area cut 
by a plane passing normal to the stirrup. 

(p. 44), Forcing Function Data Block: To input a static 
load which is generated by a forcing function in 
prescribed load increments, the following data are 
required: joint no. (col. 1-5), direction indicator 
(col. 11), the reference function number (col. 16-19), 

k (col. 20), 1.0 (SF) (col. 21-30), and LA = 0, TA = 0, 
TP = 0 (col. 31-60). In a static analysis, time is 


used as a dummy variable for defining the domain of 
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the forcing function. The sign (+ or -) of the force 
should be associated with the scale factor value 


input in the Forcina Function Data Block (p. 44). 


Modifications: 


a. 


(p. 17-18), Control Data Block: Card No. 2, Columns 21-25, 
have been changed from Al, 4X to 2Al, 3X. Column 22 will 
contain the item MPRINT = secondary print level flag. 
M = suppress paging 
S, d* = standard paging 

Card No. 3, Columns 6-10,have been changed from 
5X to 3X, I2 with the addition of the variable IINITD, 
read initial guess of displacements (activated if this 
value is negative). This option may be selected if the 
user wishes to input an initial guess of the displaced 
configuration of the structure rather than have SINGER 
make the initial guess. The initial quess displacements 
should be placed at the end of the Forcing Function Data 
Block, after the zeros card, in 8E10.0 format. They 
should be in the order of the joint numbers and the 
order of X, Y, and rotation at each joint. The joint 
displacements should be followed by the internal node 
displacements for all the elements in the order of the 
element numbers. This gives a total of initial displace- 
ments equal to the number of degrees of freedom at the 
joints plus the number of elements. The displacements 


corresponding to the joints should be normalized by the 
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average element length. Columns 61-65 have been changed to 
I5 to add the variable INKS, output print option. This 
gives the user the choice to specify the time intervals 

at which the output will be printed. For example, if 

the user wants the results printed at every 10 time 
increments, the value 10 should be placed in columns 64-65. 
This will give results according to IPRINT, which is also 
specified at every 10 time increments. The intermediate 
time steps will display the time, function value, and 
displacements in a condensed nondimensional form. 

The new subroutines DEFO and STRN, written to improve 

the directional properties of the element model (see 
section V, Volume I, and also section IV, Volume II), 
redefine the element deformation components. Since the 
existing subroutines LEAF and FAIL utilize the old 
component definitions, both should be studied for 

possible modifications before all four subroutines are 
used within the same program. 

Also, these two new subroutines do not distinguish 
between finite and infinitesimal joint rotations. The 
same transformations are used regardless of the input 
parameter ILIN (p. 17). 

Due to the function used to compute the arc 
tangent value in the new DEFO subroutine, values of 
joint rotations are valid within the range of (-1<@<r). 


Errors will be generated in this evaluation for 


rotations outside this range. 


ERROR CONTROLS 

In addition to inevitable roundoff errors, the dynamic response 
predictions of SINGER are subject to the following errors: 

1. Temporal discretization error 

2. Iteration error 

3. Spatial discretization error 

4. Quadrature error 

The temporal discretization error is caused by the representation 
of the motion during a time step in terms of a finite power series in 
time, and the iteration error is caused by the termination of the 


search process after a prescribed accuracy has been achieved. These 


two errors are controlled automatically by SINGER [2]. 


The spatial discretization error is caused by the finite element 
representation of the continuum model of the beam-column, and the 
quadrature error is caused by the numerical integration of the internal- 
energy density over the volume of the finite element. These two errors 
must be controlled by the user to assure accurate energy predictions. 

A detailed description of error controls in energy evaluations is 


presented in section IV of Volume I of this report. 


MINIMIZATION AND CONVERGENCE PROBLEMS 

Listed below are the basic components of a general problem which 
may significantly influence the solution process: 

1. Control parameters 

2. Model characteristics 


Load increments 


Saree t 
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1. Control parameters: The control parameters include the 
minimization convergence tolerance and the error measures. The state 
of convergence of the minimization process is controlled by the 
minimization tolerance parameter (SERR), a user input quantity which 
defines a limit to the stepsize used in a linear minimization. Also, 
an error value is computed which is intended to provide a measure of 
the quality of the converged state, referred to as the iteration 
error [2]. This error measurement is derived from the largest of the 
unbalanced generalized forces at the joints. In the program, this 
value is printed as the FULL TIME ERROR. For a static problem, this 
is the only error measure and it does not control the response in 
any way. For a dynamic problem, a half time error is computed, which 
is an error measure at the middle of the time step. It attempts to 
measure the truncation error due to approximate representation of the 
displacement function over a time step [2]. The size of the time 
step is controlled by the half time error. 

The convergence tolerance parameter chosen too large can produce 
inaccurate results because the process terminates before a well defined 
minimum has been reached. A lower limit is placed on the parameter 
determined by the computer word size [1]. A range of values found 


to produce acceptable results for most problems is 1. x 107 to 


Te ® 1077. It should also be noted that since the previous converged 
state is the basis for initialization of the next increment, there is 
some interaction between successive increments. However, the exact 


significance of this interaction with respect to solution accuracy 


is not well defined. 


The error measures are computed during the solution process. It 
is possible to reach a converged state with respect to the convergence 
tolerance specified but still have a poor error measure. This may be 
due to a complex energy surface in the neighborhood of the minimum 


state chosen. Experience has shown that accurate results (measured 


cat OR ee ae 
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by symmetry in the element response) are achieved when the full time 
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error is smaller than 1. x 1072; a value of 1. x 10 is questionable; 


and values larger than 1. x 107! usually produce inaccurate response 
predictions. Good results are associated with the smaller values, 


4 or smaller. 


such as 1. x 10° 
The half time error measure in dynamic problems is compared with 
a tolerance 1. x 1073 (CRITU) for determining an acceptable converged 
state. (The magnitude of this error is usually larger than the full 
time error.) The time step is reduced to a magnitude of 0.60 of the 
previous value if this convergence limit is not achieved. The minimum 
time step defined in the code is 1. x 10. It is also possible that 
the time step could be automatically increased to protect against 
a roundoff error in computations. This occurs if the half time error is 
4 less than a fixed value, 1. x 1079 (CRITL). The next time interval is 
made 1.85 times the previous interval. (For a discussion of the 
dynamic error control logic, see Ref. [2], pages 50-53.) 
2. Model characteristics: The most significant modeling 
characteristic which may influence the convergence process is the 
element distribution. The distribution chosen must accurately represent 


the continuum being modeled as reflected in how well the energy 


distribution can be reproduced. Extremely nonlinear material characteristics, 
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such as an elastic-plastic bilinear function, introduce complexities 
into the energy distribution, frequently occurring as localized energy 
concentrations at concentrated loads and totally restrained joints. 
Any region which is likely to develop large strain gradients should 
be subdivided to more accurately represent the energy concentration. 
It appears that members should be subdivided so that no one element 
contributes more than 50 percent to the energy function. The effect 
of poor element distribution can be reflected in the following ways: 
Large strain discontinuities may occur at the joints (see section VI, 
Volume I); and the accuracy of the continuum response prediction may 
be unsatisfactory. 

3. Load increments: The magnitude of toad increments can cause 
various convergence problems. A relatively large load increment may 
produce an equilibrium state so far removed from the previous state 
that convergence may be difficult to achieve, particularly if 
complexities due to nonlinear effects exist in the total energy surface. 
Also, the structure may reach a state of reduced stiffness so that a 
smal? change in load produces a relatively large change in the 
corresponding displacements. This has been observed in the steel 
wide flange beams and reinforced concrete beams loaded to failure 
(see section VI, Volume I), and also in the study of the elastica B]. 
The effect of this difficulty can usually be reduced or eliminated 
by decreasing the load increments. 

Another source of convergence difficulty occurs as the predicted 


response approaches the limit point of the equilibrium path. If the 


load is incremented to a value very near the limit point, an extremely 


large displacement state may result with a poor error measure indicated. 


power . 7 ny yen Fee 


If the load level is incremented above the limit point, the solution 


search produces such extremely large values that the computer control 


system terminates the problem. There is no automatic termination 

: check made in the program to detect this condition. 

For most structures, the choice of load increments is not critical. 

| The minimization algorithm appears to be a stable process and convergence 
is usually achieved for the load increments specified. The solution time 
may be longer for a problem containing a few large load steps as 
compared to one with several small load steps because the number of 
iterations for convergence for each large increment may be significantly 
increased. However, the solution accuracy may also be increased for 
some problems. 

Two types of normal convergence may occur. The usual and most 
desirable type of convergence occurs when the stepsize used in a linear 
minimization is reduced to the specified tolerance. The other type 
occurs when a zero stepsize has occurred in the direction of a linear 
é | minimization. This type of convergence is indicated by the printed 
message NO MOTION IN THE LINEAR MINIMIZATION. This condition may occur 
during an apparently normal minimization process. It may also occur 
when an element yields. When yielding occurs, the minimization at 
that load level is repeated with the status of the yielded member 
changed. Since yielding may have occurred only at the outer fibers 
of the element, the strains at the integration points may still be 
within the elastic range. If this is the case, no change occurs in 


the solution vector and no motion is indicated. 
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Probably the most important consideration in obtaining accurate 


predictions of nonlinear structural response is the user's understanding 
of the system behavior. This is an indispensable aid in deciding how 

to subdivide the members and simplify the mass discretizations. It 

may be necessary to obtain satisfactory results through a trial and 
error procedure. Also, the user can study a test member to determine 
the subdivision necessary to accurately represent the expected load- 


deformation states (see section IV, Volume I). 


STRAIN DISCONTINUITY 

If one element is required to represent large changes in strain 
within its length (of the order of 50 times the yield strain or 
larger), a strain discontinuity may occur with respect to an adjacent 
element which is not so severely strained. Large strain gradients 
occur in regions of highly localized deformations, such as the region 
near a concentrated load laterally applied to a reinforced concrete 
beam after the steel reinforcement begins to yield. The discontinuities 
may appear not only as differences in magnitude, but also in sign. 
The solution process does not appear to be adversely affected by these 
discontinuities. However, since the stress resultants depend on the 
stress values computed across the sections at the ends of an element, 
the strain discontinuities produce stress resultant discontinuities. 

The strain discontinuity problem can be improved by subdividing 
the member into smaller elements within the localized region. In 


this manner, each element is required to represent a smaller strain 


change along its length. Discontinuities may still occur, but their 
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magnitudes should be reduced since the energy distribution can be 


better represented. A more complete discussion with a demonstration 


of these effects is presented in section VI, Volume I, of this report. 
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SECTION III 
ADDITIONAL PROGRAM DETAILS 


Some program details have not been adequately explained in previous 
documentation. Such detailed descriptions are important for understanding 
the program structure, and they provide the basis for making additional 
modifications to the code. Program details discussed in this section 
include a description of selected variables, the organization of the 
array storage scheme called DATA and KDATA, and the procedure to be used 


for modifying the fixed storage allocation built into the program. 


DESCRIPTIONS OF SELEUTED VARIABLES 

The variables described here are primarily those which determine 
the location of Gauss points and end points for both steel wide flange 
and reinforced concrete elements. The objective is to relate the names 
used in the code to the physical location of the points in the elements. 
The Gauss points define the locations at which energy density values are 
computed, and the end points are used in part as stress points from 
which to compute the element stress resultants. These point locations 
are shown in figure 6 for a reinforced concrete element and in figure 7 
for a wide flange element (pages 31 and 34, respectively). 

The location of a point for any element type is defined by an x 


coordinate and a y coordinate. The following list shows the variable 


names associated with the locations for both element types. 


Variable Names 


X Location Y Location 
Gauss Point End Point 
Steel wide flange XPI Y@P YFIBR 
Reinforcing bar XPI YBAR YBAR 
Concrete (confined and 
unconfined) XPI Y@P YFIBR 


All elements use the same x location variable XPI as an array with 


the following form: XPI(J,M) J = section number 
M = element number 


The section number refers to the cross sections located at the three 
Gauss coordinates in the longitudinal direction (for the Gauss point 
locations), and at the ends of the element for the end point locations. 
The y location names are made identical where possible. The Gauss 
point locations for each type are named YGP, used as an array with the 


following form: YGP(1,M) I = y location number 
M = element number 


The locations on a cross section are shown in figures 2 and 3 for the 
two element types. Points defined by I = 1, 2 and 6, 7 are at the Gauss 
locations for the flange or cover thickness. The remaining three 
points defined by I = 3, 4, and 5 are at the Gauss locations for the web 
or confined concrete depth. 

The y locations for the end points use the same array name YFIBR 


y location number 


with the following form: YFIBR(1I,M) I 
M = element number 


The locations are also shown on the cross sections of figures 2 and 3. 


There is a total of 11 points for each element type. The points 
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GAUSS POINT 


LOCATIONS 


YGP(1,M) 
YGP(2,M) 


YGP(3,M) 


YGP(4,M) 


YGP(5,M) 


YGP(6,M) 


YGP(7,M) 


END POINT 


LOCATIONS 


YFIBR(1,M) 


YFIBR(2,M) 


YFIBR(3,M) 


YFIBR(4,M) 


YFIBR(5,M) 


YFIBR(6,M) 


YFIBR(7,M) 


YFIBR(8,M) 


YFIBR(9,M) 


YFIBR(10,M) 


YFIBR(11,M) 


FIG. 2 GAUSS POINT AND END POINT LOCATIONS FOR 
A REINFORCED CONCRETE SECTION. 
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GAUSS POINT END POINT 
LOCATIONS LOCATIONS 


YFIBR(1,M) 
YGP(1,M) 


YGP(2,M) 
YFIBR(2,M) 
YFIBR(3,M) 
YGP(3,M) 
YFIBR(4,M) 


YFIBR(5,M) 


YGP(4,M) YFIBR(6,M) 


YFIBR(7,M) 

YFIBR(8,M) 
YGP(5,M) 

YFIBR(9,M) 


YFIBR(10,M) 


YFIBR(11,M) 


FIG. 3 GAUSS POINT AND END POINT LOCATIONS 
FOR A WIDE FLANGE SECTION. 
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defined by I = 1, 2 and 10, 11 are the limits of the flange or cover 
thickness. The points defined by I = 3, 4, 5, 6, 7, 8, and 9 are 
determined by dividing the web or confined concrete depth by 8 and 
spacing the points evenly at the eighth points. 

The y locations for the longitudinal reinforcing bar groups use 
the array name YBAR with the following form: 


YBAR(I,M) I 
M 


bar group number 
element number 


The YBAR values are the y locations of each longitudinal reinforcing bar 
group input for an element. Figure 6a shows the points on a reinforcing 
bar with the length the same as the element length. This form is valid 
for each bar group, regardless of the input Jength, since each one is 
redefined as an equivalent group over the element length for energy 
computations. 

The determination of whether or not a reinforced concrete element 
has a confined core is based on the value stored in the array NTIES(M), 
where M is the element number. NTIES stores the number of tie (stirrup) 
groups input for a reinforced concrete element. If NTIES = U, no ties 
exist for the element and there is no confined core. The total section 
is then treated as unconfined for energy computations. If NTIES is a 


positive integer, then there are ties within the element and the energy 


is computed for both the confined core and the unconfined outer shell. 


ORGANIZATION OF DATA/KDATA STORAGE ARRAYS 


Selected data values are stored in locations organized as two 
related arrays; one stores certain floating point data (DATA), while 
the other stores related integer data (KDATA). The DATA array includes 
the following values: 

1. Curvature data 

2. Reference forcing function data 

3 Joint forcing function parameters 
4. Plastic strain data 


5. Stress history data 


The KDATA array includes the following integer values: 
1. Reference forcing function time points 
Joint forcing function parameters 


Plastic strain data location indices | 


> Ww YY 


Stress history location indices 

The following discussion provides the organizational details for these 

two arrays. 
1. Organization of the DATA Array: The DATA array contains . 

floating point data which are necessary for the solution process. It | 

contains curvature data that are used in subroutine SEEK; it includes 


function tables and forcing function data; this array also contains 


element stress-strain histories if an element has yielded (MSTAT = 3), 
or the user has input the option for storing the stress-strain historie: 


(MSTAT = 3) rather than using the changeable status (MSTAT = 2). 


General storage in the DATA array is illustrated in figure 4. 


STORAGE : 


CURVATURE DATA 


FUNCTION TABLES 


JOINT FORCING 
FUNCTION DATA 


STRESS-STRAIN 
HISTORIES FOR 
ELEMENT "1" OR "M" 


STRESS-STRAIN 
HISTORIES FOR 


ELEMENT "M" 


STRESS-STRAIN 
HISTORIES FOR 


ELEMENT "NM" OR "M" 


SUBSCRIPTS: 


LTAB | 
LTAB + 1 


LFF 
LFF?) 


LFF + 4*NFF = LP 
) 


KDATA (LPI + 1) = LP +1 


KDATA(LPI + M) 


KDATA(LPI + NM) 


FIG. 4 GENERAL STORAGE IN THE DATA ARRAY. 
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Allocation of space in the first portion of the DATA array for 
curvature data is made in subroutine FTAB. The total allocation is 
set by LTAB = NDF(NDF + 1)/2 where NDF is the number of degrees of 
freedom of the system. 

Subroutine FTAB stores the Function Table Data in the order of 
reference forcing functions by time-force points; i.e., 


DATA(LTAB + 1) = time point 1 for function no. 1 


DATA(LTAB + 2) = force point 1 for function no. 1 


DATA(LTAB + 3) = time point 2 for function no. 1 
which is continued through 
DATA(LTAB + 2*J) = force point J for function no. 1, and function 
table number 2 would begin with: 
DATA(LTAB + 2*N + 1) = time point 1 for function no. 2. 
This format continues for all function tables through the last value: 
DATA (LFF), where LFF = LTAB + 2*J(1) + ... +2*J(NTAB), where J is 
the number of time-force pairs for function table i and NTAB is the 
total number of function tables. 


Joint forcing function data are stored into DATA array by subroutine 


JFOR. The following format is used for each joint forcing function: 


DATA(LFF + 1) = A = scaling factor 

DATA(LFF + 2) = B = load addition constant 
DATA(LFF + 3) = C = time addition constant 
DATA(LFF + 4) = D = time period for sinusoidal or 


cosinusoidal functions. 
These four values are stored in this manner through the last point: 
DATA(LP), where LP = LFF + 4*NFF, and NFF is the total number of joint 


forcing functions input by the Forcing Function Data Block. 
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Element stress-strain data are stored in the array only if an 
element has yielded or the user has forced storage during input by 
setting IACT to N on the element parameter card in the Element Data 
Block. In most problems, not all elements will have yielded before 
collapse takes place. For this reason, only the elements requiring 
storage of their stress-strain data need the storage allocation, and 
these data are stored in the DATA array as each element yields. (The 
entry points for the stress and strain data are stored in the KDATA 
array.) A simple five element example where elements have yielded 
in the order of 3, 2, and 5 is shown in Figure 5. Reinforced concrete 
stress-strain data are stored in the following general form for each 
member: 

REINFORCING STEEL STRAINS: 5*NGRP(M) points 
UNCONFINED CONCRETE STRAINS: 31 points 

CONFINED CONCRETE STRAINS: 27 points 

REINFORCING STRESS HISTORIES: 40*NGRP(M) points 
UNCONFINED CONCRETE STRESS HISTORIES: 248 points 
CONFINED CONCRETE STRESS HISTORIES: 216 points 

Temporary storage is made from the DATA array in the S array 
in subroutines COEN, STEN, and OUTS. Strain and stress data are 
updated through the S array, which is contained in the common block 
FIBER, in subroutines CRET and REIN. The S array contains nine values 
which reference one point within the element. The values S(1) through 
§(8) correspond to the stress history at a particular point, and S(9) 
contains the strain at that point within the element. These values are 


directly transferred between the S and DATA arrays when MSTAT is set 
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KDATA 


SUBSCRIPT STORAGE 


LP. + T + Ro 


3 


Le? 1+ Ry 


DATA 


ELEMENT | STORAGE | SUBSCRIPT 
: Re LP +1 
S 
3 
2 2 
S 
2 
Rs 
35 


LP + 7+ R3 


LP, +1 


mers STRAIN DATA 
Sn = STRESS DATA 
LP, = LAST POINT STORED 
m = ELEMENT "m" 


NOTE: ELEMENTS HAVE YIELDED IN THE 
ORDER 3, 2, 5; ELEMENTS 1 & 4 
HAVE NOT YIELDED. 


FIG. 5 RELATION BETWEEN DATA AND KDATA; 
5-ELEMENT EXAMPLE. 
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to 3 for the element. If MSTAT is not set to 3 for the element, 
there is no transfer between these arrays at any of the points within 
the element. 

Storage in the DATA array is made in the order of Gauss points 
(integration points) and force points (end points) for each of the six 
general groups of strain and stress listed above for each element. 
These point locations are defined in the preceding subsection, DESCRIPTIONS 
OF SELECTED VARIABLES. 

Storage of the strain points of element M is made in the following 
way: 

LP = KDATA(LPI + M) - 1 is the reference point as shown in 
Figures 4 and 5. 
An example for two steel reinforcing bars is illustrated in Figure 
6a. The storage for an element containing two bars is made as follows: 
DATA (LP + 1) 
LP + 2 
LP + 3 
LP + 4 
bet 5 
LP + 6 
where 1 through 6 refer to the Gauss points (x) in Figure 6a. 
Storage continues with: 
DATA (LP + 7) 
Le +s 
LP +9 


LP + 10 


a. LONGITUDINAL STEEL BARS 


a 


b. UNCONFINED CONCRETE 


c. CONFINED CONCRETE 


FIG. 6 STRESS-STRAIN STORAGE FOR REINFORCED 
CONCRETE ELEMENTS. 
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where 7 through 10 refer to the force points (o) in Figure 6a. 

The concrete strain points are referenced by KRAIN = KDATA(LPI + M) + 
5*NGRP(M) - 1, where NGRP(M) is the number of groups of reinforcing 
bars in element M. The subscripts for the unconfined and confined 
concrete in the DATA array continue from DATA(KRAIN + 1) to DATA(KRAIN + 58) 
as illustrated in Figures 6b and 6c. 

The reinforcing steel stress histories are referenced by KRESS = 
KDATA(LPSI + M) - 1 which is illustrated in Figure 5. Eight values are 
stored for each of the Gauss points and force points in Figure 6a. In 
this example with two bars, the subscripts would continue with KRESS + 2 
through KRESS + 48 for the Gauss points (1-6) and KRESS + 49 through 
KRESS + 60 for the force points (7-10). For each reinforcing bar group 
there are 24 storage locations for Gauss points and 16 locations for end 
force points which totals 40 points for each group. 

The concrete stress histories are referenced by KRESS = KRESS + 40* 
NGRP(M). Eight values are stored at each of the points shown in Figures 
6b and 6c. Storage is made for the unconfined concrete at the Gauss 
points with subscripts KRESS + 1 through KRESS + 168 which are 21*8 
values. The force points are stored with subscripts KRESS + 169 through 
KRESS + 248 which are 5*8 + 5*8 values. Stress histories at the Gauss 
points of the confined concrete begin with the subscript KRESS + 249 and 
continue through KRESS + 320 which are 9*8 values. The force point stress 
histories continue from KRESS + 321 through KRESS + 464 which are 9*8 + 
9*8 values. 

Following is a summary of the storage in the DATA array of a 


reinforced concrete element, M, where NG = NGRP(M). 
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Strains 


Reinforcing steels: 


Gass points 


Force points 


Unconfined concrete: 


Gauss points 


Force points 


Confined concrete: 


Gauss points 


Force points 


Stress Histories 


Reinforcing steels: 


Gauss points 


Force points 


Unconfined concrete: 


Gauss points 


Force points 


Confined concrete: 


Gauss points 


Force points 


Values Stored 


3*NG 
2*NG 


18 


Values Stored 


24*NG 
16*NG 


168 
80 


72 
144 


Total 


3*NG 
5*NG 


5*NG + 21 
5*NG + 31 


5*NG + 40 
5*NG + 58 
Total 


24*NG 
40*NG 


40*NG + 168 
40*NG + 248 


40*NG + 320 
40*NG + 464 


Storage in the DATA array for a steel wide flange element is 


as shown in Figure 7. 


from KRAIN + 22 through KRAIN + 43. 
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done in a similar way as that for a reinforced concrete element. 
indexing for the strain data begins with KRAIN = KDATA(LPI + M) - 1. 
The Gauss point (x) subscripts range from KRAIN + 1 through KRAIN + 21 
The subscripts for the force point strains range 


The stress histories are referenced 


en 


FIG. 7 STRESS-STRAIN STORAGE FOR WIDE 
FLANGE ELEMENTS. 
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by KRESS = KRAIN + 43. The subscripts for the Gauss point stress his- 
tories range from KRESS + 1 through KRESS + 168 which are 21*8 values. 
The stress history subscripts for the force points range from KRESS + 169 
through KRESS + 344 which are 11*8 + 11*8 values. 


{ 
| The following is a summary of storage in the DATA array of a 


wide flange element (Figure 7): 
| Strains Values Stored Total 
Gauss points 21 21 
Force points 22 43 
Stress Histories Values Stored Total 
Gauss points 168 168 
Force points 176 344 


2. Organization of the KDATA Array: The KDATA array contains 
integer data necessary for the solution process and index locations for 
the stress-strain histories in the DATA array. It is first stored by 
subroutine FTAB with function table data. The variable LTABI, which 
is the load table index, is the locator for each function table input 

into the program through the Function Table Data Block. This variable 
) is initialized to zero in ABET. The function table data are stored in 
the following format: 


KDATA(LTABI + 1) 


LTABI + ] 

KDATA(LTABI + 2) = number of time points for each 
reference forcing function through 
the last forcing function table: 
KDATA(LFFI), 

where LFFI = LTABI + 2*NTAB 


and NTAB is the total number of reference forcing functions. 
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Forcing function data are stored in KDATA in subroutine JFOR through 
input in the Forcing Function Data Block. Storage is in the following 
format: 

KDATA(LFFI + 1) = loaded joint number 
KDATA(LFFI + 2) = reference joint number 
KDATA(LFFI + 3) = direction 
KDATA(LFFI + 4) = reference function number 
KDATA(LFFI + 5) = last forward time point used for 
each joint forcing function through 
the last input forcing function: 
KDATA(LPI), 
where LPI = LFFI + 5*NFF and NFF is the total number of joint forcing 
functions used. 

The index or subscript values that define the locations of the 
element stress-strain histories in the DATA array are set in subroutine 
STOR. The following format is used: 

KDATA(LPI + 1) = start of strain data for element 1 

KDATA(LPI + 2) = start of strain data for element 2, and 
for each element, where KDATA(LPSI) is the 
starting point for the strain data of element 


NM and LPSI = LPI + NM. 
KDATA(LPSI + 1) 


start of stress history data for element 1 


KDATA(LMAXI ) start of stress history data for element NM 


where LMAXI = LPSI + NM = LPI + 2*NM. 


MODIFICATION OF FIXED STORAGE ALLOCATION’ 

When a problem is sufficiently large or computer storage is critical, 
various arrays within SINGER must be modified to accommodate these 
restrictions. The key to successful storage modification is an under- 


standing of the several key variables which describe in general the 
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problem size in question. These variables are defined as follows: 


Variable Definition Current Value 
NDFD Maximum number of degrees of freedom 90 
NJD Maximum number of joints 50 
NLD Maximum number of leaf springs 20 
NMATD Maximum number of materials 9 
NMAX Maximum index of DATA array 10000 
NMAXI Maximum index of KDATA array 500 
NMD Maximum number of members 45 


Every change of storage within SINGER can be accomplished by changing 
those arrays which depend on these seven variables. Other program 
variations are possible, but a complete reprogramming would be necessary. 
These include having greater than 10 groups of reinforcing steel in one 
member, varying the Gaussian quadrature system in the program, and varying 
the number of elemental distortion components (varying the number of 
internal nodes). 

Table 1 gives a list of the common storage data arrays and the 
variable arrays that must be changed. Table 2 gives other variables that 
must be dimensioned differently. Table 3 shows where the different common 
blocks are found within the subroutines. Table 4 gives example sizes 
for linear and nonlinear problems of the KDATA and DATA arrays. 

Given a fixed storage requirement, the largest problem that can be 
solved by SINGER would consist of one where the material response is 
strictly linear in nature. All arrays, even though they are not used 


by a particular problem, should have a minimum value of |. 


TABLE 1 - COMMON VARIABLE ARRAYS 


NECESSARY 
VALUE 


CURRENT 
VALUE 


VARTABLE 
ARRAY 


LOCATION 


DATA 10000 NMAX 


COMMON 


KDATA 500 NMAXI 


COMMON/CONBK COAREA 4,45 4 ,NMD 
5 ,NMD 


SIGMA 


5,45 


IP 45 NMD 
45 
1Q 45 NMD 


NLD 


COMMON/ELEMENT 


NLD 


IQL 45 


MATR NMD 


45 


MATW 45 NMD 


MBAR 10 ,NMD 


NMD 


MCODE 


MSHEAR NMD 


NMD 


MSTAT 


NMD 


MTIES 


MTYPE NMD 


NGRP NMD 


NMD 


NSPACE 


NMD 


NTIES 


COMMON/FIBER DENS NMATD 


EC 


NMATD 


oe oa © 


NMATD 


NMATD 


NMATD 


n 
oat 
”n 


COMMON/JOINTS ACC 


DER 

RESENG 

IDFI 

IDFII 

COMMON/MEMBER AGRP 
ATIES 


NMATD 
NMATD 

8 ,NMATD 
8 ,NMATD 
-8,NMATD 
NMATD 
NMATD 
NMATD 

3 NID 

3 NID 

3 NUD 
3,NUD 

3 NUD 

3 NID 

3 NUD 

3 NUD 
3,NJD 
3,NJD 
NJD 
3,NJD 
NJD 

3 NUD 
3,NUD 
NDFD 
NDFD 

10 ,NMD 


6 ,NMD 


4 we 


= hee 


a 


PDP 
SPRING 
STIES 

TFWF 
TWWF 
UDM 
URM 
XBEG 
XBEGM 
XBEGS 
XL 


45 
45 
10,45 
45 

45 

45 

45 

45 
10,45 
45 

45 

45 

45 
7,45 
5,20 
7,45 
45 

45 

45 

45 
10,45 
45 
6,45 
45 
5,45 
10,45 
7,45 


NMD 
NMD 
10,NMD 
NMD 
NMD 
NMD 
NMD 
NMD 

10 ,NMD 
NMD 
NMD 
NMD 
NMD 

7 ,NMD 
5 ,NLD 
7,NMD 
NMD 
NMD 
NMD 
NMD 

10 ,NMD 
NMD 

6 ,NMD 
NMD 

5 ,NMD 
10 ,NMD 


7 ,NMD 


Bist Sa 


COMMON/SAVEBK 


COMMON/SEEKBK 


YFIBR 
YLDS 
XDM 
PDF 

DISM 


SAVACC . 


SAVAXL 
SAVCRY 
SAVMOM 
SAVSHR 
SAVSRP 
SAVSRO 
SAVXDJ 
SAVVEL 
SVSTRN 
SVSTRS 
DEFOR 
STPSIZ 
GRAD 
GRADI 
DELTAG 
DIRECT 
DIAG 


11,45 
45 

45 
7,45 
45 
3,50 
2,45 
2,45 
2,45 
2,45 
3,20 
3320 
3,59 
3,50 
12,45 
12,45 
90 

90 

90 

90 

90 

90 

90 


11,NMD 
NMD 
NMD 

7 ,NMD 
NMD 

3 ,NJD 
2,NMD 
2 ,NMD 
2 ,NMD 
2 ,NMD 
3,NLD 
3,NLD 
3 ,NJD 
3 ,NJD 
12 ,NMD 
12 ,NMD 
NDFD 
NDFD 
NDFD 
NDFD 
NDFD 
NDFD 


NDFD 
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TABLE 2 - OTHER VARIABLES 


SUBROUTINE VARIABLE (CURRENT SIZE) 

ABET SOLN(90) 
NDFD=90 

NJD=50 

NLD=20 

NMATD=9 

NMAX=10000 

NMAXI=500 

NMD=45 

ADYN SOLN(90) 
ASAN SOLN(90) 
CUTS KIND(45) 
DELT SOLN(90) 
ELIN KIND(45) 
ERRS SOLN(90) 
JFOR TERROR(50,5) 
MATD TERROR(9,16) 
POTE SOLN(90) 
REJO IRESTR(3,50) 
TERROR(5 ,50) 

JNUM(50) 

SEEK DEFORM(90) 
CURV(4095) 


NECESSARY VALUE 


NDFD 
NDFD 
NJD 
NLD 
NMATD 
NMAX 
NMAXI 
NMD 
NDFD 
NDFD 
NMD 
NDFD 
NMD 
NDFD 
NJD,5 
NMATD , 16 
NDFD 
3,NJD 
5 ,NJD 
NJD 
NDFD 


NDFD* (NDFD+1 ) 
—— 


TABLE 3 - SUBROUTINES AND COMMON BLOCKS 


COMMON BLOCKS TO CHANGE 


SUBROUTINE 


COMMON/ ELEMENT 
COMMON/ FIBER 
COMMON/JOINTS 
COMMON/MEMBER 
COMMON/SAVEBK 
COMMON/SEEKBK 


4 
co 
= 
oO 
oO 
io 
— 
oOo 
= 
= 
oO 
oO 


¥841939S/NOWWOD 


GIAVS/NOWNOD 


43 aW3W/NOWWOD 


SLNIOL/NOWWOD 


44g14/NOWWOD 


LN4W314/NOWWOD 


€NO9/NUWNOD 


NOWWO9 


SUBROUTINE fi 


8433S /NOWNOD ae | 
adAVS /NOWNOD hy ol 
Yd aWaW/NOWNWOD a ee 
SLNIOC/NOWNOD mele sy 


Y4a14/NOWWOI >< oe 
~< ~< i 


~< 


4N00/NOWNOI 
NOWWO9 =< =< 


SUBROUTINE 
BLOCK DATA 


i RC ei me 


TABLE 4 - SIZE OF KDATA AND DATA ARRAYS 


LINEAR PROBLEM 
ARRAY NECESSARY SIZE CURRENT VALUE 
KDATA LPI = (2*NTAB)+(5*NFF)+(2*NM) 500 
DATA LP = (NDF*(NDF+1)/2) NE 249(2)+4*NFF 10000 

NDF = Number of dagrene 06 freedom 


yon 


NFF = Number of joint forcing functions 
NM = Number of members 


NTAB = Number of forcing function tables 

J(i) = Number of time-force pairs 
NONLINEAR PROBLEM 
ARRAY NECESSARY SIZE CURRENT VALUE 
KDATA LPI 500 
DATA LP+NMC*[40*NG+464 ]+NMWF*344 10000 


NG = Number of groups of reinforcing 


steel in member 
NMC = Number of R/C members 
NMWF = Number of WF members 


TYPICAL PROBLEMS 
Problem A 


10 members 


40 degrees of freedom maximum 


2 forcing functions 


Problem B 


50 members 
200 degrees of freedom 
10 forcing functions 


6 time-force pairs for each function 6 time-force pains for each 


8 joints loaded 


All R/C members with three groups 


of reinforcing steel 


KDATA 
Linear Problem SIZE 
Problem 64 
Problem B 220 


Nonlinear Problem 
roblem 64 
Problem B 220 


20 joints loaded 
All R/C members with three 
groups of reinforcing steel 


DATA 
SIZE 


20300 


6716 
49500 


SECTION IV 
NEW SUBROUTINE DOCUMENTATION 


Three new subroutines were developed and one was significantly 
altered for efficiency and improved accuracy. Subroutine COEN was 
modified by removing a segment of computations which could be more 
efficiently done in a separate subroutine during the input sequence. 
The new subroutine, CONC, is described in this section as well as the 
revised organization of COEN. The description of subroutine COEN is 
available in the Program Document of a previous contract report [4]. 
In addition, two new subroutines were written to replace DEFO and 
STRN as required by the details of the new finite element model (for 
development details, see section V, Volume I, of this report). The 
description details for these two subroutines are basically the same 
as the original documentation [4]; only the computational details 
have been changed. The source listings for both DEFO and STRN are 
included in this section since they are currently introduced into the 
program at the option of the user, and they are not included in the 


complete program listing in the Appendix of this volume. 
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Name of Subroutine - CONC(M) 


Description of Subroutine 


a. 


Purpose and Uses. This subroutine determines the cross- 
sectional area of the confined concrete and the areas 

of the sides, top, and bottom of the unconfined concrete 
for each element. The integration constants of the 
Gauss-Legendre quadrature are incorporated into these 
area computations for use in subroutine COEN. This 
routine also alters points 5 and 6 of the confined 
concrete stress-strain curve when ties or stirrups 

are within the element. 

Description of Input. The element number (M) is input 
through the argument list. Element properties and 
stress-strain data are obtained from COMMON data 
storage. 

Description of Expected Output. The output quantities 
are the factored cross-sectional areas of concrete for 
each element (M) and the altered points of the stress- 
strain curves. 

Limitations/Restrictions. The cross-sectional areas 


are factored for a Gauss-Legendre quadrature. 


Relationship to Other Routines. This routine is called 
by subroutine ACIN. 


START 


COMPUTE AREAS OF 
CONCRETE SECTIONS 


REMOVE CONCRETE AREA 
DISPLACED BY STEEL 


ALTER POINTS 5&6 
OF STRESS-STRAIN CURVE 


COMPUTE SLOPES 4&5 
RETURN 


SUBROUTINE COEN 


START 
RETRIEVE INDEXES FOR 
STRESS-STRAIN HISTORIES 
CALCULATE STRAINS AT 
INTEGRATION POINTS 
<> 


YES 


: 
| 


a 


ESTABLISH CONFINED CONCRETE 
STRESS-STRAIN DATA 


DETERMINE INTERNAL ENERGY 
DENSITY AT EACH POINT 
CALL CRET 


SUM CONTRIBUTIONS OF 
ALL CONFINED POINTS 
UPDATE MEMBER 
3 END STATE DATA 
a DETERMINE INTERNAL ENERGY DENSITY AT 
q EACH POINT OF UNCONFINED CONCRETE 
CALL CRET 
SUM CONTRIBUTIONS OF 
ALL UNCONFINED POINTS 
UPDATE MEMBER 
END_STATE DATA 


ee 


10 
20 


40 
30 


SUBROUTINE DEFO(M) 

COMMON/ELEMET/ICARD,IP(45), IPL(20),1Q(45),1QL(20),1DUM1(1125) 
COMMON/JOINTS/DUM2(1500) ,X(50) ,XDJ(3,50),¥, (50) ,0UM3(480) 
COMMON/LEADBK/AVDM,AVGL ,DUM4 (5) ,1D1(20) ,DT EPS, 1D2(20) ,DUM5(11) 
COMMON/MAINBK/IDUM2(42) ,NPRT, IDUM3(6) 
COMMON/STRNBK/SRP(4) ,SRQ(4) ,UX,UY,UZ,XLEN,AREA,ZZ1, IMAT 
IF(M.GT.0) GO TO 10 

iLs = TABS(M) 


= IPL(ILS) 
= IQL(ILS) 
LEN = 1.£0 
GO TO 20 
I = IP(M) 
J = 1Q(M) 
DX1 = X(J) - X(1) 
DX2 = Y¥(J) - Y(I) 
DUT = (XDJ(1,J)-XDd(1,1))*AVGL 
DU2 = (XDJ(2,J)- XDJ(2, 1) )*AVGL 
DX1P = DX1 + DUI 
DX2P = DX2 + DU2 


DXLEN = SORT(DX1**2+DX2**2) 
DXLENP = SORT (DX1P**2+DX2P**2) 
ADX = OXT*OX2P - DX2*DXiP 
BDX = DX1*DX1P + DX2*DX2P 
IF(BOX.NE.0.E0) GO TO 30 
WRITE (NPRT,40) 
FORMAT (///10X,43H****** ARGUMENT BDX (DEFO) IS ZERO ******//) 
RETURN 
RBANG = ATAN2(ADX,BDX) 
sr = 0.5E0* (DXLENP-DXLEN) 

= XDJ(3,I) - RBANG 
UZ = Wg J(3,J) - RBANG 


SUBROUTINE STRN(M,X,Y,STRAIN) 
COMMON/MEMBER/ DUMMY 1 (5365) ,XDM(45) ,DUMMY2 (360) 
COMMON/ STRNBK/SRP(4) ,SRQ(4) ,UX ,UY,UZ,XLEN,AREA,ZZ1, IMAT 
ETA = Y/XLEN 

= 2.E0/XLEN*X - 1.E0 
ALPHA = (UY**2 - 0.5EO*UY*UZ + UZ**2)/15.E£0 


BETA = (UZ**2 - UY**2)/16.E0 

STRAIN = (2.E0*UX/XLEN+ALPHA) + (BETA-4.E0*XDM(M) )*T 
] - ETA*((3.E0*T-1.E0)*UY + (3.E0*T+1.£0)*UZ) 
RETURN 


END 
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APPENDIX 


PROGRAM SOURCE LISTING 


The complete listing of the SINGER code is included in this 
section with the exception of the two new subroutines DEFO and STRN 
developed to implement the logic of the new finite element model i 
(section V, Volume I). The listing for these two subroutines is 


included in section IV of this volume. 
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